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ABSTRACT 

The search for substellar companions around stars with different masses along the main sequence is critical to understand the different 
r- H processes leading to the formation of low-mass stars, brown dwarfs, and planets. In particular, the existence of a large population of 

low-mass stars and brown dwarfs physically bound to early-type main-sequence stars could imply that the massive planets recently 
imaged at wide separations (10-100 au) around A-type stars are disc-born objects in the low-mass tail of the binary distribution and 
are thus formed via gravitational instability rather than by core accretion. Our aim is to characterize the environment of early-type 
main-sequence stars by detecting substellar companions between 10 and 500 au. The sample stars are also surveyed with radial 
velocimetry, providing a way to determine the impact of the imaged companions on the presence of planets at £ 10 au. High contrast 
and high angular resolution near-infrared images of a sample of 38 southern A- and F-type stars have been obtained between 2005 
and 2009 with the instruments NaCo on the Very Large Telescope and PUEO on the Canada-France-Hawai'i Telescope. Direct and 
saturated imaging were used in the J to K s bands to probe the faint circumstellar environments with contrasts of ~ 5 x 10~ 2 to 10" 4 
at separations of 072 and 1", respectively. Using coronagraphic imaging, we achieved contrasts between 10~ 5 and 10~ 6 at separations 
> 5". Multi-epoch observations were performed to discriminate comoving companions from background contaminants. This survey 
is sensitive to companions of A and F stars in the brown dwarf to low-mass star mass regime. About 41 companion candidates 
were imaged around 23 stars. Follow-up observations for 83% of these stars allowed us to identify a large number of background 
contaminants. We report the detection of 7 low-mass stars with masses between 0.1 and 0.8 M in 6 multiple systems: the discovery 
of a M2 companion around the A5V star HD 14943 and the detection of HD 41742B around the F4V star HD 41742 in a quadruple 
system; we resolve the known companion of the F6.5V star HD 49095 as a short-period binary system composed by 2 M/L dwarfs. 
We also resolve the companions to the astrometric binaries i Crt (F6.5V) and 26 Oph (F3V), and identify a M3/M4 companion to the 
F4V star o Gru, associated with a X-ray source. The global multiplicity fraction measured in our sample of A and F stars is > 16%. 
jj^ This has a probable impact on the radial velocity measurements performed on the sample stars. 
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1. Introduction sive enough for triggering the runaway accretion of hydrogen 

(see, e.g., Mordasini et al. 2009). On the other hand, ~ 3% of 

Stars do not form alone, as they are often found in multiple sys- lanets ^ detected fe direct ■ { searches at j g 

H terns, commonly consisting of two or more stellar components, tions ^ 1Q ^ ef erentially Mund early-type stars. Most of 

brown dwarfs, or planets. The various properties of these sys- these { A lanets seem tQ haye ft diflferent Qri ^ name , ^ 

O terns in terms of mass ratios and separations call for different fr tation of the prot0 pi ane tary gaseous circumstellar disc 

formation mechanisms. At the lower end of the mass distnbu- ( Dodson-Robinson et al. 2009; Kratter, Murray-Clay & 

tion, more than 400 extrasolar planets have been detected, most Youdin 2my ^^ there couW bg tWQ iWe hs of langt 

> f T ( ~u 8 ° %) vel ° clmetrlc ™ys and in most cases formatiori5 dep ending on the stellar mass and planet separation. 

• rH the planet-host stars are of solar type. This technique brought 

/\ observational evidence that giant gaseous planets at separations The core-accretion theory (Pollack et al. 1996) is currently 

?-j S 6 au, which correspond to the largest revolution periods prefered to the gravitational instability theory (Boss 1997) as 

probed by radial velocimetry today, are preferentially forming the origin mechanism of most planets detected by velocimetry 

through rapid accretion of gas on pre-existing rocky cores, mas- because the latter mechanism is not expected to produce giant 

planets at distances ^ 10 au (Rafikov 2005). Meanwhile, obser- 

Send offprint requests to: D. Ehrenreich vational support has been gathered in favour of the core accretion 

* Based on observations made with ESO Telescopes at the Paranal theory: the correlation between the planet frequency and the host 

Observatory under programme IDs 076.C-0270, 081.C-0653, and sta r metallicity (Santos, Israelian & Mayor 2001, 2004) and the 

083.C-0151, and on observations obtained at the Canada-France- emerging fact that planets with an intermediate mass between 

Hawai l Telescope (CFHT) which is operated by the National Research KT . , jo»> j-i i»u j-^j 

_ ., r „ j , T . „ T ■ i j n • j mi ■ r Neptune s and Saturn s are exceedingly rare are both predicted 

Council ol Canada, the Institut National des Sciences de 1 Univers of ... , . , r ™ r , . . 

the Centre National de la Recherche Scientifique of France, and the ln thls theoretical frame. The recent extent of velocimetnc sur- 

University of Hawai 'i ve y s t0 low-mass main sequence stars (M dwarfs; Mayor et al. 

** The fun version of this paper with the appendices in available on line at 2009a,b) yielded new key evidence in favour of this formation 

http://www-iaog.obs.ujf-grenobie.fr/~dehrenre/articies/afsurvey/ mechanism: the fact, for instance, that low-mass planets (about 
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the mass of Neptune) are frequent while gas giants are seldom 
found around low-mass stars observed by radial velocity or mi- 
crolensing techniques (Bonfils et al. 2007; Sumi et al. 2010). On 
the contrary, the frequency of giant planets should not be im- 
pacted by the stellar mass in the frame of gravitational instabil- 
ity, as long as circumstellar discs are massive enough to become 
unstable (Boss 2006). 

In fact, the total mass of the circumstellar disc where planets 
form, is supposed to scale with the mass of the central star. An 
increase of the disc total mass can lead to an enhanced growth 
rate for protoplanetary cores in the disc midplane (Ida & Lin 
2004), leading to a larger number of massive planets around 
early-type stars than around Sun-like stars. Hence, there should 
be an observational correlation between planet occurrence and 
stellar mass (Laughlin, Bodenheimer & Adams 2004; Ida & Lin 
2005; Kennedy & Kenyon 2007). Surveying early-type, main- 
sequence stars hotter and more massive than the Sun should al- 
low to test such a correlation. 

Radial velocimetry can be used to detect planets at short 
periods around main-sequence early-type stars, as shown by 
Lagrange et al. (2009a). These authors used the HARPS spec- 
trograph installed at the ESO 3.6-m telescope in La Silla (Chile) 
to survey a sample of 185 stars. They measured on each star the 
achievable detection limits, taking into account the 'jitter' ' level 
of each object. These authors estimated that planets with periods 
up to 100 days could be found around ~ 50% of the surveyed 
stars. Constraining the presence of planets at larger separations 
however requires a different detection technique. 

In this respect, direct imaging is a powerful tool used in com- 
plement to radial velocimetry for detecting companions at sep- 
arations typically <; 10 au. It is sensitive to intrinsically bright 
objects: low-mass stars, brown dwarfs, and young giant plan- 
ets, and bring essential insights on the possible origins of these 
objects (Chauvin et al. 2010). Direct near-infrared imaging al- 
lowed observers to detect several kinds of substellar companions 
such as brown dwarfs (e.g., Nakajima et al. 1995; Lowrance et 
al. 1999, 2000; Chauvin et al. 2005a), objects at the transition 
between brown dwarfs and planets (e.g., Chauvin et al. 2005b), 
or planetary-mass objects such as 2M1207b in orbit around a 
brown dwarf (Chauvin et al. 2004, 2005c); all of which are likely 
not formed through core accretion like planets detected by ve- 
locimetry, but rather as multiple stars via a cloud fragmenta- 
tion process or as brown dwarfs through disc fragmentation (see, 
e.g., Lodato, Delgado-Donate & Clarke 2005). These detections 
nevertheless show that star-forming mechanisms can be efficient 
down to planetary masses. 

Recent breakthroughs in high-contrast imaging enriched this 
picture, as giant planets were directly imaged around Fomalhaut 
(Kalas et al. 2008), HR 8799 (Marois et al. 2008), and/? Pictoris 
(Lagrange et al. 2009b, 2009c, 2010); all these stars being young 
and early main-sequence A stars. The discovery of a compan- 
ion to fi Pic, a 9 ± 3-Mj planet with a semi-major axis of 
~ 8 to 15 au, implies that giant planets can indeed form in 
~ 10 Myr. According to theoretical prediction of Kennedy & 
Kenyon (2007), this particular planet is close enough to its host 
star and could have formed in-situ by core accretion. This is not 
so clear for Fomalhaut b, which is located 115 au away from 
its star, meaning that a particular migration scenario is required 
if the planet formed closer to the star via core accretion (Crida, 
Masset & Morbidelli 2009), or for the three planets or brown 
dwarfs located at 68, 38, and 24 au from HR 8799. Gravitational 



1 The excess of scatter in radial velocity measurements resulting from 
inhomogeneities on the stellar surface and stellar envelope pulsations. 



instability, rather than core accretion, seems more suited to 
explain the existence of such massive planets on wide orbits 
(Dodson-Robinson et al. 2009). The picture is however not that 
simple, as the gravitational instabilities leading to the fragmen- 
tation of massive circumstellar discs seem to produce objects 
with masses typically above the deuterium-burning planetary- 
mass limit. In fact, according to Kratter, Murray-Clay & Youdin 
(2010), atypical disc conditions are required to form planetary- 
mass object via this mechanism. These authors suggest that if 
these planets formed this way, they must lie in the low-mass tail 
of the disc-born binary distribution. In this case, a larger number 
of brown dwarfs or low-mass stars (such as M stars) should be 
found around A-type stars at distances of 50-150 au. 

In this frame, we perfomed a deep imaging survey of early- 
type A and F stars included in the velocimetric survey of 
Lagrange et al. (2009a). This study would allow us to test 
whether low-mass stars and brown dwarfs commonly cohabit 
with massive main-sequence stars, bringing new constraints on 
the origin of the massive planets imaged around these stars. In 
addition, this imaging survey could help determining the impact 
of stellar multiplicity on the presence of closer-in planets, de- 
tected in parallel with radial velocimetry. Binarity is indeed an- 
other critical parameter for the theory of formation and evolution 
of planets. In particular, the separation of the binaries could im- 
pact the way giant planets form, either by core accretion or disc 
instability (Zucker & Mazeh 2002; Eggenberger, Udry & Mayor 
2004; Desidera & Barbieri 2007; Duchene et al. 2010), and the 
models of binary discs now predict observational effects. Mayer 
et al. (2005), for instance, predict a lack of planets around each 
component of binary systems with separations % 100 au because 
of instabilities in the circumstellar discs. The subsequent dynam- 
ical evolution also depends on the properties of the star and the 
outer companion (see, e.g., Rivera & Lissauer 2000). 

Only adaptive optics (AO) deep imaging allows to test prop- 
erly the presence of massive substellar companions. A num- 
ber of studies set to probe the existence and the impact of 
such objects to exoplanetary systems detected by velocimetry 
has been previously undertook (Patience et al. 2002; Luhman 
& Jayawardhana 2002; Chauvin et al. 2006, 2010; Mugrauer, 
Seifahrt & Neuhauser 2007; Eggenberger et al. 2007). For in- 
stance, Eggenberger et al. (2007) searched for bright long-period 
companions around 130 G- and K-type stars and measured a bi- 
nary fraction of (8.8 ± 3.5)% in a sub-sample of 57 stars with 
planets detected with radial velocimetry at separations of 0"8- 
6"5 (i.e., 30-250 au, whereas they found a higher binary fraction 
of (12.3±3.2)%inacontrol sub-sample of 73 stars with no planet 
found. Previous works were also targeted on particular objects, 
such as Vega or £ Virginis (Hinkley et al. 2010 and references 
therein), in search of substellar companions. 

In this work, we are focusing on the close environment of a 
sample of austral A- and F-type stars, which we mainly observed 
from the southern hemisphere using the AO system NaCo in- 
stalled at the ESO Very Large Telescope (VLT) on Cerro Paranal 
(Chile). Some austral stars were also observed from the north- 
ern hemisphere using PUEO, the AO bonnette of the Canada- 
France-Hawai 'i Telescope (CFHT) on Mauna Kea (USA). In a 
few cases, we also used NaCo to observe stars with positive de- 
clinations. This paper reports on these observations, which are 
described in Sect. 2. The data reduction is detailed in Sect. 3, we 
present our analysis in Sect. 4 and discuss the results in Sect. 5. 
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2. Star sample and observations 

2.1. Sample selection and biases 

The present survey is the imaging part of the radial-velocity sur- 
vey described in Lagrange et al. (2009a) that was designed to 
detect planets around early-type stars with the HARPS spectro- 
graph at the ESO 3.6-m telescope in La Silla (Chile). The sur- 
vey is limited to dwarfs with spectral types ranging from F7 to 
B8 2 . Figure 1 presents the stars in a Hertzsprung-Russell dia- 
gram while Fig. 2 shows the number of stars observed per spec- 
tral type. In total, we surveyed 38 stars, including 16 F-, 19 A-, 
and 3 late B-type stars. 

The velocimetric and imaging surveys are also volume- 
limited, with distance limits set at 33 and 67 pc for the F0-F7 
and B8-A9 dwarfs, respectively. Figure 3 gives the number of 
stars observed per 5-pc distance bin. The difference in distance 
limits would allow us to have roughly the same number of A 
and F stars in the sample. Note, however, that one early F star 
(HD 4293) with a distance of 66.6 pc is included in the survey. 

Finally, the few stars present in our sample that are not 
included in Lagrange et al.'s (2009a) sample are part of the 
northern radial-velocity survey in progress with the SOPHIE 
spectrograph (Bouchy et al. 2009) at the 1.93-m telescope in 
Observatoire de Haute-Provence (France). The properties of all 
stars in the sample can be found in Table 1. 

While more details about the sample selection are given in 
Lagrange et al. (2009a), we emphasize three biases, (i) First, 
spectroscopic binaries and close visual binaries with separations 
smaller than 5" known at the beginning of the velocimetric sur- 
vey were excluded from the target list. Hence, these were not 
considered in the imaging survey as well, (ii) In addition, the 
present imaging survey is in fact strongly biased towards 'in- 
teresting' targets, i.e., stars around which companion candidates 
(CC) were detected during early epochs. This strategy was set 
mainly in order to compensate for the limited amount of observ- 
ing time devoted to the imaging programme, (iii) Finally, poor 
atmospheric conditions and technical problems that were experi- 
enced especially during observing runs VI and VII, prevented us 
from obtaining second-epoch observations for 13% of those tar- 
gets with companion candidate(s) identified during a first epoch; 
atmospheric conditions making the AO correction loop unstable 
also prevented us from obtaining the highest achievable contrast 
- through the use of the coronagraph - for 10% of the targets. 

2.2. Observations 

Data were recorded during seven observing runs (or epochs) per- 
formed between January 2005 and August 2009, at the VLT and 
CFHT, where we totalized 7 + 8 nights of observations, respec- 
tively. The observing dates, instruments, and filters used during 
each epoch are summarized in Table 2, while the detailed ob- 
serving set-ups for each target are given in Table 3. 

2.2.1. VLT/NaCo 

We used the NAOS-CONICA instrument (NaCo) set on the VLT 
Unit Telescope 4 (Yepun) to benefit from both the high im- 
age quality provided by the Nasmyth Adaptive Optics System 
(NAOS; Rousset et al. 2003) at infrared wavelengths and 
the good dynamics offered by the Near-Infrared Imager and 
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Fig. 1. Colours of stars in the imaging sample. Stars with spectral types 
F, A, and B are represented red-, green-, and blue-filled circles, respec- 
tively. 
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2 Velocimetric detection limits for stars earlier than B8 does not fall 
into the planet domain. 



Spectral type 
Fig. 2. Spectral types of the stars in the sample. 



Spectrograph (CONICA; Lenzen et al. 2003) detector, in or- 
der to study the close circumstellar environment of 38 early- 
type stars. The NaCo Shack-Hartmann visible wavefront sensor 
was chosen to perform the AO corrections on these bright tar- 
gets, used as self-references. NaCo allowed us to perform direct 
imaging as well as coronagraphic imaging in order to improve 
the image contrast. We used the coronagraphic mode consist- 
ing in a Lyot stop in the pupil plan of the telescope, combined 
with a occulting mask of diameter © = 077 inserted in the focal 
plane. According to the atmospheric conditions, we used dif- 
ferent broad- and narrow-band filters whose properties (central 
wavelength A c , bandpass AA, and transmission T) are listed in 
Table 4. For instance, poor atmospheric conditions degrade the 
Strehl ratio S ; using a long-wavelength filter (typically K s ) al- 
lows to compensate for this degradation since the value of S 
also scales with wavelength as S — exp-(2nu)/A) 2 , where u> is 
the root-mean-square deviation of the wavefront. 

Two objectives were employed in order to optimize the point 
spread function (PSF) sampling in these different bandpasses. 
The S13 camera has a field of view of 14 x 14 arcsec 2 and a 
mean plate scale of 13.21 mas pixeL 1 ; it was preferentially used 
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Table 2. Epoch dates and astrometric calibrations for all observing runs. 



Epoch 


UT Date 


Nights 


Instrument/Camera 


Filter 


Plate scale 
(mas) 


True north 
O 


I 


2005-01-27 


4 


PUEO/KIR 


Bry 


35.06 ±0.1 


0±0.1 


11 


2005-11-06 


3 


NaCo/S27 


J 


27.06 ± 0.06 


0.02 ±0.1 








NaCo/S27 


K, 


27.06 ± 0.06 


0.02 ±0.1 


111 


2007-01-28 


2 


PUEO/KIR 


Bry 


35.08 ±0.03 


-2.39 ± 0.04 








PUEO/KIR 


K 


34.87 ± 0.03 


-2.44 ± 0.03 


IV 


2007-11-16 


2 


PUEO/KIR 


Bry 


35.04 ±0.03 


-2.63 ± 0.02 








PUEO/KIR 


Fen 


34.83 ± 0.01 


-2.64 ± 0.02 


V 


2008-08-20 


2 


NaCo/S27 


K s 


27.08 ± 0.03 


0.04 ± 0.09 








NaCo/S13 


K s 


13.27 ±0.04 


-0.04 ± 0.22 








NaCo/S13 


H 


13.26 ± 0.04 


-0.02 ±0.21 


VI 


2009-04-26 


1 


NaCo/S27 


K, 


27.04 ± 0.02 


-0.41 ± 0.06 








NaCo/S13 


H 


13.19 ±0.06 


-0.88 ± 0.37 








NaCo/S13 


J 


13.19 ±0.05 


-0.74 ± 0.26 


VII 


2009-08-27 


1 


NaCo/S27 


K, 


27.06 ± 0.02 


-0.36 ± 0.05 
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Fig. 3. Number of stars per 5-pc distance bin. The F, A, and B stars are 
represented by the histograms filled with tight diagonal red stripes, di- 
agonal green stripes, and horizontal blue stripes, respectively. The dot- 
ted line stands for the whole star sample. 



when observing in J and H bands. When observing in the K s 
band, the S27 camera was chosen: it offers a 28 x 28 arcsec 2 field 
of view and a mean plate scale of 27.06 mas pixel -1 . For each 
epoch, precise plate scales were redetermined using astrometric 
calibrators (see Table 2). 

Our observing strategy consists, for each target, in obtain- 
ing a first-epoch image and, when companion candidates are de- 
tected, to perform a second-epoch observation in order to test 
their status: either field objects (background contaminants) or 
comoving objects (physically bound to the targeted star). The 
discrimination between these two possibilities is allowed by the 
stellar proper and parallactic motions in the plane of the sky. This 
is illustrated in Fig. 4. Hence, depending on motion amplitudes, 
a few months to a few years are necessary between the first- and 
second-epoch observations. For instance, for two observations 
performed three years apart, and a precision on the measured po- 
sitions of individual point sources on each epoch of ~ 10-20 mas 
(see Table 5), we can roughly estimate that a 3-cr discrimina- 
tion between comoving and non-comoving sources is possible 
for typical minimum proper motions of ~ 15-30 mas yr~' . 

A target is first observed both in direct imaging and corona- 
graphic modes. For epochs V, VI, and VII, the CoroObsAstro 



observing template was used: this template consists in acquiring 
first some coronagraphic exposures on the object without any 
neutral density filter (Full_Uszd setting, where the pupil area is 
simply reduced by the presence of the Lyot stop, which blocks 
14% of the light). Then, a neutral density filter (ND_Short) is in- 
serted in the focal plane while removing both the coronagraphic 
mask and the Lyot stop in the focal and pupil plan, respectively. 
A direct image is then taken, the object being at the same exact 
location as behind the mask. This allows for retrieving the star 
position behind the mask, enabling precise astrometric measure- 
ments of separation with companion candidates only visible in 
coronagraphic mode. The star is then jittered (by steps of ~ 5") 
on 4 different locations on the detector in order to obtain the 
background sky estimation for the direct images. The last offset 
sends the telescope ~ 30" away from the target star, and the oc- 
culting mask is used again - the neutral density filter is removed 
once more - as 5 jittered exposures are taken on the sky in order 
to estimate the sky background for coronagraphic images. This 
template was not available during our first NaCo observing run 
(epoch II), hence the precise location of the star behind the mask 
is not known, and the astrometry is therefore less precise than 
for later observations. In all cases, the same amount of time was 
spent on the object and on the sky. 

A systematic astrometric shift is induced when switching the 
neutral density filter back and forth: we used pinhole calibrations 
recorded at each epoch to ensure that this shift is small: typically 
+0.1 pixel and ±0.01 pixel along the detector x- and y-axes, re- 
spectively. This shift is included in our error budget within the 
<j s term of Eq. (3). 

When a companion candidate can be detected without the 
coronagraph, the star is reobserved at a later epoch, only in direct 
imaging mode. NaCo observing templates ImgGenericOf f set 
and ImgAuto Jitter were used for this purpose. 

2.2.2. CFHT/PUEO 

PUEO (Rigaut et al. 1998) is the AO bonnette mounted on the 
3.6-m CFHT. It was used in combination with the near-infrared 
camera KIR (Doyon et al. 1998) which has a field of view of 
35 X 35 arcsec 2 , and a mean plate scale of 34.8 mas pixel -1 . 
Depending on atmospheric conditions and source brightnesses, 
we used broad- and narrow-band filters listed in Table 4. 

We performed unsaturated direct imaging to investigate as 
close as possible to the star, while the detector was saturated in 
order to increase the contrast at larger separations (there is no oc- 
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Table 4. Properties of the filters used during the survey. 



Filter 


A c (jim) A/1 (jum) 


T(%) 


VLT/NaCo 


J 


1.27 


0.25 


78 


H 


1.66 


0.33 


77 


NB212 


2.122 


0.022 


55 


NB 2 . 17 


2.166 


0.023 


52 


K s 


2.18 


0.35 


70 


CFHT/PUEO 


H 


1.632 


0.296 


85 


Fen 


1.650 


0.015 


70 


K' 


2.115 


0.350 


92 


H 2 (v = 1-0) 


2.117 


0.020 


65 


Bry 


2.169 


0.020 


70 


K 


2.198 


0.336 


90 



3.2. Coronagraphic data 

The reduction procedure makes use of the eclipse utility li- 
brary (ESO 1996-2002; Devillard 1997). For each observing 
template on a target, coronagraphic exposures are sorted by eye. 
Good exposures (i.e., exposures where the star is well occulted 
and the AO correction loop is closed) are then concatenated into 
a data cube. A coronagraphic sky is created as the median of all 
dithered coronagraphic exposures on the sky, and is corrected 
twice for bad pixels detected (i) on the reduced flat field and (ii) 
on the median sky once the first bad-pixel rejection as been ap- 
plied. The cleaned median sky is then subtracted to each plane 
of the data cube, which are next divided by the flat field. The 2 
bad-pixel maps are again used on the flat-fielded data cube. The 
final reduced image is extracted as the median of the data cube. 



culting mask on PUEO). A special care was taken to ensure that 
the unsaturated exposures could be used as references to mea- 
sure the relative position and photometry between the star and 
the companion candidate(s). Images are jittered on the detec- 
tor so that the sky background can be retrieved. The observing 
method is otherwise similar to that described above for NaCo. 
The 6 stars from our sample with PUEO epochs have all been 
observed with NaCo, and are thus included in this study for con- 
sistency. 



3. Data reduction 

3.1. Calibrations 

Flat fields have been taken on the twilight sky for each observ- 
ing run at VLT (usually within a few days from the observing 
nights). At CFHT, dome flat fields were taken for every filter at 
least once per run. 

A well-known field in the Orion Trapezium centered on 
the star 8 l Orionis C was used for astrometric calibrations 
(McCaughrean & Stauffer 1994). One observation per fil- 
ter/camera set-up used was taken at each observing epoch to 
determine the mean plate scale and the true north orientation. 
On each observation, the exposure time was adjusted to satu- 
rate the 2 or 3 brightest stars of the field so that from 15 to 40 
stars (depending on the size of the field, from 13" with the S13 
camera of NaCo to 35" with PUEO) have a sufficient signal-to- 
noise ratio (> 5). The position of the centroid (s x , s y ) of each 
non-saturated star on the reduced image of the field is measured 
in pixels; these values are then compared to the position on the 
sky {p x ,Py) in arcseconds reported by McCaughrean & Stauffer 
(1994) with the following equations: 



p x - x Q + p x s x cos 0„ + PySy sin 0„ 
Py = yo ~ PxS x sin 6„ + PyS y sin 6 n 



(1) 
(2) 



where p x (p y ) is the plate scale along the x- (y-) axis, 6 n is the 
orientation of the detector on sky, and xq and yo are offsets giv- 
ing the correspondance between the absolute positions (it is only 
used to solve the equation). A Levenberg-Marquardt minimiza- 
tion is then applied to find the plate scale solution. The errors on 
the plate scale and orientation are given by the standard devia- 
tion of the corresponding parameter. The mean plate scales and 
orientations are reported in Table 2. 



3.3. Direct imaging data 

The idea is similar to the reduction of data taken in corona- 
graphic mode, except that all the steps are performed by the 
jitter routine, which realigns dithered direct images via a 
cross-correlation algorithm and does all the cosmetic tasks. An 
additional (but similar) reduction is devoted to the first direct im- 
age taken after the coronagraph is removed (see Sect. 2.2.1); this 
reduced image will be used as an astrometric reference for the 
position of the star behind the mask. 



4. Analysis 

4.1. Extraction of astrometric and photometric parameters of 
candidate companions 

Candidate companions are spotted by a close eye-examination of 
each star, and their absolute position on the detector determined 
to a ~ 1 -pixel precision. The position of the star is similarly 
estimated in unsaturated direct images; it is estimated from the 
astrometric direct images for coronagraphic observations. These 
positions are then used as inputs by a dedicated IDL wrapper pro- 
gramme, that successively performs three astrometric and pho- 
tometric estimations. An aperture estimation (i), where the aper- 
ture is centered on the input position, returns the centroid po- 
sition as the barycentre of the light within the aperture, and an 
arbitrary magnitude estimation from the integration of all pix- 
els within the aperture, weighted by the flux of the reference 
PSF within the same aperture. The position (Ax, Ay) and magni- 
tude Am of each companion candidate are then given relatively 
to those of the star. A Gaussian fitting (ii) using the IDL rou- 
tine GCNTRD is also applied to each input position and returns 
relative separations. Finally, (iii) a deconvolution algorithm us- 
ing a Levenberg-Marquardt minimization of the log-likelyhood 
function (Veran 1997) is employed: it basically consists in a PSF 
fitting within the pupil plane. Procedures (i) and (ii) work well 
when the objects are well separated on the image, and give sim- 
ilar results than procedure (iii). However, they do not give accu- 
rate results for tight systems (typically tighter than the aperture 
radius used). In these cases, our estimations rely on procedure 
(iii). 

The astrometric precision is ultimately limited by the astro- 
metric calibrator, and the precision cr p (in mas) obtained on the 
estimation of the plate scale p (in mas pixeU 1 ) at each epoch 
(which are listed in Table 2). Yet, our procedure for extracting 
the astrometric parameters of an image is also characterized by 
an uncertainty <r s on the measured separation s (both in pixels). 
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The total uncertainty on the measured separation p (in mas) is 
thus 



(3) 



The values of cr s depend on various parameters such as the actual 
separation of the components as well as their contrast. In the fol- 
lowing, we choose to take <x s - 0.5 pixel as a conservative value. 
In fact, this choice does not impact our conclusions concern- 
ing moderately separated systems. Meanwhile, this value can be 
thought of as representative since it is clear that the real precision 
we can achieve does not allow us to resolve the most difficult 
cases, such as the very tight companion to HD 43940. 

Since the stellar flux is always estimated from the direct (un- 
saturated) image, both the direct and coronagraphic images are 
normalized by their respective exposure times when the pho- 
tometry is performed for objects present in NaCo coronagraphic 
images. The magnitude difference is then increased by 4.6 mag 
to take the transmission difference introduced by the use of a 
neutral density filter into account. In the case of NaCo corona- 
graphic imaging, the use of a slightly undersized aperture due to 
the presence of the Lyot stop in the pupil plan leads to an ad- 
ditional uncertainty of ~ 4% on the fluxes measured when the 
occulting mask is on, with respect to the full aperture available 
when it is off (Boccaletti et al. 2007). 

As for the separation, the determination of the contrast be- 
tween two components is also impacted by their respective 
brightnesses and separation. We can actually estimate our pho- 
tometric precision from the set of dispersion of the measure- 
ments taken on the same objects at different epochs: the median 
dispersion obtained is 0.4 mag (45% on the measured fluxes). 
Astrometric and photometric measurements for all companion 
candidates imaged are given in Table 5. Stars without any de- 
tected companion candidates are listed in Table 6. 

4.2. Companionship 

The companionship, i.e., the fact that a candidate companion 
is comoving with the primary star in the vincinity of which it 
has been imaged - which is indicative of a physical associa- 
tion bound by gravity - is appreciated by measuring the evo- 
lution of the separation p = \[{Ax 2 + Ay 2 ) and the position angle 
9 - arctan Ax/ Ay (from north to east) between the candidate and 
the star over several epochs of observation. The values of Ax and 
Ay are corrected for the orientation 9„ of the detector north with 
respect to the true north (given in Table 2 at each epoch). The 
shift in the position (right ascension a and declination <5) of a 
companion candidate between two epochs i and j is 



Acvj^j - pj sin 0j - pi sin #, 
A8i^,j - pj cos 9j - pi cos 9 t 



(4) 
(5) 



These measured shifts are compared to the theoretical appar- 
ent motion the object should have on the sky, given the stellar 
proper motions listed in Table 1, the parallactic motion, and as- 
suming, as standardly done, that the candidate is a motionless 
background contaminent. Figure 4 shows two examples of com- 
panionship determination. The positions of all companion can- 
didates obtained during several observing runs are represented 
in Fig. Al (see the Appendices 3 ). 



The full version of this paper with the appendices in available on line at 
http : //www-laog . obs . u j f-grenoble . f r/~dehrenre/articles/af survey/ 



For observations performed at N epochs, it is possible to cal- 
culate the probability Pbkg that a companion candidate is a back- 
ground object, using a x 2 probability test of 2N - 2 degrees of 
freedom: 
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Xbkg 
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(6) 



where Aa* 



and AS* 



are the theoretical shifts in a and 6 in 



the position of the star between epochs 1 and i, taking the stellar 
proper motion and the parallactic motion into account. 

The probability of observing a value of x 2 that is larger than 
that obtained with Eq. (6) for a random sample of N observa- 
tions with v = 2N — 2 degrees of freedom is the integral of the 
probability density of a x 2 -distribution (Bevington & Robinson 
2003), 
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(7) 



where the T(n) function is equivalent to the factorial function 
n\ extended to nonintegral arguments, and (x 2 ) represents the 
possible values of x 1 within the integral sum. The probability 
Pbkg 1S obtained from Eq. (7); practically, we are using the IDL 
chisqr_pdf function. A status of background contaminant is 
assigned to each object for which Pbkg > 0.01%. This probabil- 
ity is given for each object in Fig. Al. This first test allows to 
reject 13 out of 41 detected point sources (32%), identified as 
background contaminants (labelled '(B)ackground' in Table 5). 
Because of systematic effects in the position measurements, 
due for instance to variations in the stellar proper motions, dis- 
crepancies between the real values and those from the literature, 
or the fact that background objects have a non-negligible proper 
motion on the sky, we calculated Pbkg < 0.01% for some objects 
that are evidently not comoving with the primary star. This is, 
for instance, the case of HD 68456 CC#1 or HD 91889 CC#1 
(see Fig. Al). For this reason, we additionally define the proba- 
bility P cmv , also given in Fig. Al, that a companion candidate is 
comoving with the star, using a test similar to Eq. (6), 



Acmv 



-s 



;=i L 
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0"Aff, 



ASi 



0~A<5, 



(8) 



The result is injected into Eq. (7). Note, however, that the value 
of -Pcmv can only be used to give a hint of a physical associa- 
tion. In fact, this probability does not take into account the or- 
bital motion of a true comoving companion. Hence, an object 
with very small values of Pbkg and P cmv (< 0.01%) is not nec- 
essarily a background contaminant because of orbital motions: 
this is typically the case of HD 49095 CC#1 and CC#2, which 
are revolving around each others while clearly following the star 
over the 4 observing epochs. In the following, we consider that 
an object with i\kg < 0.01% and P cmv > 0.01% is comoving 
with the star: 4 companion candidates out of 41 detected point 
sources (10%) are confirmed comoving objects according to this 
test (they are labelled '(C)omoving' in Table 5). Consequently, 
an object with Pbkg < 0.01% and P cmv < 0.01% could be con- 
sidered to be a background contaminant. On this basis, we can 
confidently reject 6 additional point sources (15%). 

Another tricky case is that of HD 41742 CC#1 (P bkg = 0%, 
P cmv = 0%), observed during 4 epochs: as can be seen in Fig. Al 
the second-epoch measurement is far from that expected of a co- 
moving object; it is also far from that of a background contami- 
nant. Meanwhile, measurements recorded at all other epochs are 
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Table 6. Sample stars around which no companion candidates were detected. 
* At CFHT, saturated imaging was used instead of coronagraphic imaging. 



# 


Star name 


Date 


Filter Camera 


DIT (s) 


x NDIT x NEXP 












Direct imaging 


Coronagraphic imaging 


1 


HD 2834 


2005-11-06 


K 


S27 


0.8x5x 10 


3.0x9x22 


2 


HD 3003 


2005-11-07 


K 


S27 


0.8x5x 10 


6.0x3x30 


3 


HD 4293 


2005-11-08 


K 


S27 


1.0x5x 10 


10.0 x 3 x 20 


4 


HD 10939 


2005-11-06 


K 


S27 


0.4x5x 10 


5.0 x 6 x 20 


7 


HD 16754 


2005-11-07 


K 


S27 


0.4x5x 10 


5.0 x 6 x 20 


8 


HD 19545 


2008-08-20 


K 


S27 


2.0x5x 10 




12 


HD 29992 


2008-08-20 


K 


S27 


0.8 x 13 x 10 


2.0 x 15 x 10 


13 


HD 31746 


2005-11-07 


K 


S27 


0.6x5x 10 


8.0x5x 15 


23 


HD 112934 


2007-01-28 


K 


KIR* 


1.0 x lx 18 


30.0 x 1 x 30 






2008-08-20 


K 


S27 


3.0x4x 10 




24 


HD 116568 


2008-08-20 


K 


S27 


l.Ox lOx 10 




28 


HD 186543 


2008-08-20 


K 


S27 


1.0 x 12 x 10 


5.0x6x 10 


30 


HD 200761 


2008-08-21 


K 


S27 


l.Ox 10 x 10 








2009-08-27 


K 


S27 


0.35 x 30 x 10 


l.Ox 30 x 10 


31 


HD 209819 


2008-08-21 


K 


S27 


0.35 x 30 x 10 








2009-08-27 


K 


S27 


0.35 x 30 x 10 


l.Ox 30 x 10 


34 


HD 216627 


2008-08-21 


K 


S27 


0.35 x 25 x 10 








2009-08-27 


NB 2 12 S27 


2.0x30x 10 


l.Ox 5 x 10 


37 


HD 223011 


2008-08-20 


K 


S27 


5.0x2x 10 





well compatible with a comoving object (see HD 41742 CC#2 
for comparison). This object is considered to be comoving be- 
cause the measurement made during the second epoch for all 
three companion candidates to this star seem to be outliers. If 
we reject the measurements obtained at this epoch for this star, 
then CC#1 is comoving while CC#2 and CC#3 indeed are back- 
ground contaminants. 

Finally, two point sources observed at several epochs remain 
with an ambiguous nature (labelled '(A)mbiguous' in Table 5): 
HD 16555 CC#1 which is a true 'in-between' case (see Fig. Al), 
and HD 43940 CC#1, for the reasons developped in Sect. 4.4.2. 
The 13 other imaged point sources (32%) have only been ob- 
served once; their status is set to '(U)ndefined' in Table 5. 

4.3. Detection limits 

For each image, we derive 6-<x detection limits in the form of a 
2-dimensional map. The detection limits are calculated by mea- 
suring for each pixel the noise (standard deviation) in a 5 x 5 
box centered on that pixel (see also Lagrange et al. 2009c). This 
operation is performed with the IDL routine IMAGEJVARIANCE 
(by M. Downing). The coronagraphic images are normalized to 
the exposure time of the associated direct image before their 
detection limits are calculated; the neutral density filter trans- 
mission is also taken into account. The classic 1 -dimensional 
detection limits are then radially extracted from the 2D maps: 
the ID detection limit at a separation p from the star position 
at (Aa, A<5) - (0, 0) is the azimuthal median of all pixels in an 
annulus of mean radius p with 5-pixel thickness. 

These ID detection limits can be used for determining the 
overall properties of our survey and comparing it to other sur- 
veys. However, since the contrast is mainly limited by the pres- 
ence of speckles at close separations, we believe that the 2D 
maps are more accurate than the ID limits and should rather be 
used when refering to particular objects in the survey. At large 
separations from the central stars, the contrast is rather readout- 
noise limited, and both ID or 2D limits could be used. 

Figures 5 and 6 show the 2D- and ID-detection limits 
achieved both in direct imaging and coronagraphic modes for 



two stars from our survey: an early A-type star (HD 2834; Fig. 5) 
and a late F-type star (HD 91889; Fig. 6) are chosen as typi- 
cal examples. Maps of the detection limits for all surveyed stars 
are available in Appendix B. The overall ID detection limits in 
the K s band are presented in Fig. 7. Unlike the direct imag- 
ing observations, the coronagraphic imaging mode did not use 
a neutral density filter, and therefore enabled improved sensi- 
tivities for wider-separation (greater than 2") sources. With the 
coronagraph in place, we were able to employ large Detector 
Integration Times (DITs), while still avoiding detector satura- 
tion. The larger DITs helped reduce detector read-out noise, 
which currently limits the wide-separation sensitivities. At sep- 
arations greater than ~ 5", we measured ~ 4-mag improvement 
between direct imaging and coronagraphic imaging modes. 

Since the different spectral types of stars included in the sur- 
vey have (i) different distance cut-off values (see Fig. 3) and (ii) 
different median ages, it is interesting to plot the overall detec- 
tion limit for each spectral type (A, F, and B) as the median of 
all detection limits. Note that since there are only 3 B stars in 
the survey, the median limit for B stars (blue curves) cannot be 
considered as representative as for A and F stars (red and green 
curves, respectively). 

In this scope, the performance in contrast is slightly better, 
typically by SS 0.5 mag, for F-type stars than for A-type stars. 
Using direct imaging mode, we achieve typical contrasts of 4, 6, 
and 8.5 mag at 072, 075, and > 1" from the star, respectively. 
The use of an occulting mask allows to obtain better contrast 
than with direct imaging from ~ 175. A maximum contrast of 
~ 13 mag is reached at separations larger than 5" from the star. 

4.4. Companion candidates detected 

Companion candidates are detected at all separations above the 
detection limits. Table 5 presents the observing dates, integration 
times, projected separations, position angles, and magnitude dif- 
ferences of the 41 point sources we identified throughout our ob- 
serving campaign. Based on the method presented in Sect. 4.2, 
we were able to determine that 20 objects (49% of sources de- 
tected) were background sources and 6 objects (15% of sources 
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p (arcsec) 

Fig. 7. Median 6-<x detection limits in the K s band for F- (red), A- 
(green), and B-type (blue) stars, in direct imaging (plain lines) and coro- 
nagraphic (dashed lines) modes. The coronagraphic detection limits are 
better than the direct imaging ones from 2" on, for the reasons detailed 
in Sect. 4.3. The measured separations p and magnitude differences AK S 
of companion candidates confirmed through multi-epoch observations 
and those with an undefined status are reported by filled and empty cir- 
cles, respectively. The confirmed companion candidate to i Crt, which 
was only observed in the H band, is marked by the red square assuming 
a(H - K) colour close to 0. The dash-dotted line atp = 0735 indicates 
the occulting mask radius. 
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Fig. 8. Absolute magnitudes vs. projected separations. The legend is 
the same than for Fig. 7. The absolute magnitudes M Ks and projected 
separations are calculated for each object given its apparent magnitude 
m Ks in the K s band and distance d from Earth. The 6-cr detection limits 
shown for F- (red), A- (green), and B-type (blue) stars, in direct imaging 
(plain lines) and coronagraphic (dashed lines) modes are the medians of 
the detection limits Mg s = f(p(d),d,m Ks ) calculated for all stars of a 
given spectral type. 



a selection of remarkable cases of companions with undefinied 
status and some that we rejected as comoving objects. 

4.4.1. Confirmed comoving companions 

For companions comoving with their primary stars, the projected 
physical separation in astronomical units can be simply extracted 
from the knowledge of the separation measured on the sky and 
the parallax (usually, the Hipparcos parallax; ESA 1997). The 
masses of companions can be retrieved using mass-magnitude 
relationships yielding from evolutionary models. The detection 
limits of our survey of A and F stars (Fig. 7) show that we 
are sensitive to companions in the brown-dwarf to low-mass 
star regimes. Thus, we chose the evolutionary models for solar 
metallicity low-mass stars and the associated mass-magnitude 
relationships provided by Baraffe et al. (1998) to determine the 
masses of our detected comoving companions. Pratically, the 
masses were obtained for given absolute magnitudes and ages 
by interpolations into Baraffe et al.'s tables (VizieR On-line Data 
CatalogueJ/A+A/337/403). 

HD 14943 — A companion with an absolute magnitude of 
M K — 6.6 mag has been identified at two different epochs 
273 (141 au) away from the A5V star HD 14943 (HIP 11102). 
According to Baraffe et al.'s (1998) evolutionary models for 
solar-metallicity low-mass stars, this companion should be a 
0.4 M star, assuming an age of 850 Myr, determined in the 
frame of Su et al.'s (2006) survey for debris discs around A stars 
with Spitzer/MlPS. Hence, HD 14943B, shown in Fig. 9a, could 
be a main sequence M2 star (Cox 2000). 4 

HD 41742 — This high proper-motion star is part of a physical 
ternary hierarchical system (Multiple Star Catalogue — MSC 
— , Vizier Online Data Catalogue J/A+ AS/1 24/75, Tokovinin 
1997). We identify CC1 to component B (CCDM J06046- 
4504B; see Fig. 9b), a 0.79 M Q star with a period of ~ 1 400 yr. 
The mass given in the MSC was determined from the measured 
B - V = 1.01, typical of a K0/K2 star and according to Cox 
(2000). The two other companion candidates (CC2 and CC3) are 
background objects. The component C of the system catalogued 
in the MSC is too far to enter NaCo or PUEO field-of-view 
(pc - 19677). Furthermore, the primary star was shown to be a 
spectroscopic binary according to the radial velocity curves ob- 
tained by (Lagrange et al. 2009). The imaged companion is not 
responsible for the radial velocity variations; hence, HD 41742 is 
a member of a quadruple system. We measured for HD 41742B 
(CC1) a projected separation of 579, i.e., 159 au and an abso- 
lute magnitude M K = 4.0 mag, which translates into a 0.8-M o 
star using Baraffe et al.'s (1998) mass-magnitude relationships, 
in good agreement with the mass previously determined. This es- 
timation was made assuming an age of 3.7 Gyr, according to the 
Geneva-Copenhagen Survey of Solar Neighbourhood (GCSSN, 
Vizier Online Data Catalogue V/130, Holmberg, Nordstrom & 
Andersen 2009). 



detected) were in fact comoving with their primary stars. These 
confirmed companions are shown in Fig. 9 and their physical 
properties summarized in Table 7. Because a second epoch ob- 
servation could not have been obtained for all targets, 15 point 
sources are still of undefined nature. In the following, we detail 
the properties of all confirmed companions as well as those of 



HD 49095 — This high proper-motion star with radial veloc- 
ity variations (Lagrange et al. 2009) is referenced as a binary by 
Makarov & Kaplan (2005) based on the difference of proper mo- 



4 According to Drilling & Landolt's calibration of MK spectral types 
(in Cox 2000, p. 389, Table 15.8). This star exhibits radial velocity vari- 
ations (Lagrange et al. 2009). 
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Fig. 10. The F6.5V star HD 49095 and its binary comoving companion 
composed of CC1 (blue circle) and CC2 (red circle). The axis are la- 
belled in projected separation with respect to the primary star centroid 
position (white star). The positions of objects on this image are retrieved 
using the deconvolution algorithm described in Sect. 4.1. 
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Fig. 11. Motion of the binary companion CC2 of HD 49095 (zoom from 
Fig. 10) relatively to CC1 (blue point) between 4 different epochs cov- 
ering 3.5 yr. The best-fit circular and face-on orbit is figured, with a 
semi-major axis of 2.27 au. 



tions measured between two different catalogues ('AyU binary'), 
namely, the Hipparcos and Tycho-2 catalogues. Here, we im- 
age the binary component and resolve it as a tight binary in 
which orbital motion is detected over several epochs. The sys- 
tem is shown as seen on 2005-11-07 in the K s band with NaCo 
in Figs. 9c and 10 and the revolution of the B (CC1) and C (CC2) 
components around their center of mass is shown over 3.5 years 
(4 epochs) in Fig. CI. The position of HD 49095C with respect 
to HD 49095B was retrieved using the deconvolution algorithm 
described in Sect. 4.1 and is plotted in Fig. 11. For the two 
VLT images (taken on 2005-11-07 and 2009-04-26) where the 
HD 49095BC system is well resolved, we assumed 1-cr position 
errors of 0.25 pixel. A position error of 0.5 pixel has been as- 
sumed for CFHT images (taken on 2007-01-27 and 2007-11-26) 
where the system is barely resolved. The tight system projected 
separation is seen to vary between 1.4 and 2.6 au. Yet, because 
we do not have enough points to fit a precise orbit, we simply 
calculated a best-fit circular orbit, assuming the system is seen 
face-on, and determined a semi-major axis of 2.27 au (y 2 = 3.3 
with 1 degree of freedom). HD 49095C accomplishes roughly 
60% of this circular orbit in 3.5 yr, which leads to a revolution 
period of 5.7 yr. Using Kepler's third law, we can then determine 
a total mass for the BC system of 0.1 1 M G . This is about three 
times smaller than the total mass determined photometrically us- 
ing Baraffe et al.'s (1998) models (0.175 + 0.13 = 0.305 M ) 
given the measured magnitudes in the H band of Mhb = 8.5 and 
M H c = 9.0 mag and the age of 2 Gyr provided by the GCSSN. 
We attribute this discrepancy to the fact that the orbit is not well 
constrained with 4 epochs. 

HD 101198— The i Crateris (HD 101198; CCDM J11387- 
1312A) entry in the CCDM mentions a companion (KUI 58B) at 



p = 1'.'4 and 6 = -134° in 1934, with a visual magnitude 7=11. 
The star is recorded as a probable single star in the Double 
and Multiple Systems Annex of the Hipparcos Catalogue (ESA 
1997), yet Makarov & Kaplan (2005) as well as Frankowski, 
Jancart & Jorissen (2007) see it as a A/y binary. Were KUI 58B a 
background contaminant, it would lie in 2009 atp w 13'.' 45 and 
9 « -141 deg, so that we should be able to detect it on NaCo K s 

images, unless its V-K colour index is below 1.5. The other 

possibility is that KUI 58B is indeed bound to t Crt. In this case, 
it is likely that CC1, shown in Fig. 9d, is that object. This would 
imply for CC1 colour indexes V — K » 11 - 7.25 = 3.75 and 
V - H x 1 1 - 7.44 = 3.56, typical of a M3 dwarf (Ducati et al. 
200 1) 5 . The mass derived from the measured Mh = 5.2 mag us- 
ing evolutionary model and assuming an age of 4.8 Gyr (from the 
GCSSN) is 0.57 M G , more typical of a M0 or a late K star (Cox 
2000). Furthermore, the primary shows radial velocity variations 
(Lagrange et al. 2009). 

HD 153363— 26 Ophiuchi (HD 153363) is referenced as 
a radial-velocity variable (Lagrange et al. 2009) as well as 
an astrometric binary (Proper motion derivatives of bina- 
ries, Vizier Online Data Catalogue J/ AJ/1 29/2420, Makarov & 
Kaplan 2005; HIP binaries with radial velocities, Vizier Online 
Data Catalogue J/A+A/464/377, Frankowski, Jancart & Jorissen 
2007), yet no orbit information is available in the litterature. The 
binary component expected by astrometric measurements could 
well be CC1. In fact, we determined that the motion of CC1 seen 
between August 2008 and April 2009 cannot be due to the com- 
bination of the stellar proper motion and the parallactic motion, 
which are for this low-proper motion star of similar amplitudes. 



' Here, we assume that the binary has the same age as the primary 
5 Gyr) and is therefore on the main sequence. 
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Table 7. Physical properties of comoving companions. 



# 


Star 


CC# 


Projected Separation 

(AU) 


Mass 
(Mo) 


5 


HD 14943 


1 


141 


0.40 


15 


HD 41742 


1 


159 


0.80 


17 


HD 49095 


1 


32.3 


0.175 


17 


HD 49095 


2 


31.9 


0.13 


22 


HD 101198 


1 


25.0 


0.57 


25 


HD 153363 


1 


11.3 


0.70 


36 


HD 220729 


1 


14.9 


0.30 



It can be seen in Fig. Al that the position of CC1 in April 2009, 
imaged in Fig. 9e, is not compatible with that of a background 
field object. Hence, the motion of CC1 can rather be explained 
by the orbital motion around 26 Oph. We observed 26 Oph B 
in three band passes and measured AK S =2.1, AH = 2.7, and 
A J = 3.1. Given the 2MASS JHK magnitudes of 26 Oph and 
its distance of 33.3 pc, the companion has absolute magnitudes 
M Ks - 4.2, M H - 4.9, and Mj = 5.5. These magnitudes are 
compatible, within the given uncertainty range of ±0.4 mag, 
with a 1.0-1.3-Gyr-old 0.7-M Q object according to Baraffe et 
al.'s (1998) model. Assuming the J - K s and H - K s colours 
are each overestimated by ~ 0.4 mag (a l-cr variation), they in- 
deed correspond to the colours of an early K dwarf (Ducati et 
al. 2001). Additional fainter and farther components of this sys- 
tem may exist; however, the confirmation of their companion- 
ship would require a second epoch with coronagraphic observa- 
tion. Meanwhile, we can already exclude CC2, also seen on the 
direct images, as an additional physical companion to 26 Oph. 

HD 220729— Lagrange et al. (2009a) identified o Gruis 
(HD 220729) as a spectroscopic binary and noted that it was 
associated to a ROSAT source by Suchkov, Makarov & Voges 
(2003). In this work, we confirm the binary nature of this star 
as we resolve it as a tight (074 or 14.9 au) and contrasted 
(AK = 4.7) binary system, as can be seen in Fig. 9f. We esti- 
mate the absolute magnitude of o Gru B to Mk - 6.7 for an 
age of 1.1 Gyr (GCSSN), which yields a mass of 0.3 M , i.e., 
a M3/M4 main sequence star. The companion detected in the 
adaptive-optics survey is not responsible for the radial-velocity 
shifts, hence o Gru is a triple system. 



position of the object measured in 2005 corresponds to that of 
a background object, non-comoving with the primary star. We 
obtain Pbkg = 24.9%. The lack of reported uncertainties in the 
CCDM does not allow to firmly conclude on the companionship 
status of CC1/CD-49 1541B; we consider it as a probable back- 
ground contaminant. We are nevertheless resolving it as a tight 
apparent binary and, since a second tight binary system is also 
present in the surroundings of rj l Pic, the system will need to be 
monitored in the future. 



HD 43940 — Clearly seen in 2005-November observations in 
the K s band, the close companion candidate to HD 43940 (p = 
0722) is harder to resolve from the primary star in 2008-August, 
both in K s and H bands. This may be due, however, to the 
poorer observing conditions compared to 2005-November. A fi- 
nal attempt to characterize this new tight system was done on 
2009-April, but the companion could barely be resolved from 
the primary star in the J band. HD 43940 has very low proper 
motions in the direction of CC1 (6 a = -5.51 mas yr , S s - 
-22.63 mas yr -1 ), but since the separation is small, it is hard to 
disentangle between the two possibilities: (i) the system is ap- 
parent, CC1 is a nearly motionless background contaminant ap- 
parently getting closer to HD 43940 because of this star proper 
motions, or (ii) the system is physically bound and the difficulty 
to resolve it as time passes results from the orbital motion of 
CC 1 around the primary. 

HD 158094— The B8 star S Arae (HD 158094) is part of an op- 
tical (Torres 1986) double system (CCDM 17311-6041) whose 
secondary component lies out the field of view (p = 4774, 
9 - 313°). Concerning the star itself, Spitzer/MlPS 24-^m pho- 
tometry indicates an age of 125 Myr and nonexistent infrared ex- 
cess (Rieke et al. 2005). Low-mass companions were searched 
for around this star by Hubrig et al. (2001) using the AO sys- 
tem ADONIS on the ESO 3.6-m telescope in La Silla. These 
authors report detection limits of AK ~ 5 mag and AK ~ 9 
for separations of ~ 2"-5" and ^ 4", respectively. This means 
CC3 (p = 1273587, 9 = 62.0°, AK = 8.2) could have been de- 
tected by Hubrig et al. (2001): it is slightly above their detection 
limit, and close to the edge of, yet within, their 24" x24" field of 
view, centered on the star, in 1999-March. Note that if CC3 were 
a background contaminant, then given the star proper motion, it 
would have been closer to star in 1999 (p = 1 173909). 



4.4.2. Companion candidates with undefined or ambiguous 
status 

HD 32743— 77 1 Pictoris (HD 32743, HIP 23482) is catalogued 
in the CCDM as a binary with an astrometry recorded in 1946 
of p - 1076 and 9 - -162°, and a photometry of AV = 
7.4 mag (no associated uncertainties on these parameters are 
reported). We observed this object during a single epoch on 
2005-11-07, and obtained the following astrometry for CC1: 
p = 127212 ± 070.0285 and 9 = -176.32 ± 0.12°. We mea- 
sured infrared apparent magnitudes of J — 12.6 (AJ = 7.8) 
and K s = 11.8 (AK S = 7.4) for CC1, which seem compati- 
ble with the visible apparent magnitude (V = 13) of the known 
companion to rj l Pic. Hence, it is likely that CC1 is the binary 
companion to rj l Pic and catalogued as CD-49 1541B. To es- 
timate whether CC1 is comoving with rj l Pic, we produced a 
proper motion plot similar to those of Fig. Al and that is pre- 
sented in Fig. Dl. We assumed arbitrarily large error bars on 
the 1946 astrometric parameters tabulated in the CCDM. The 



4.5. Notes on some remarkable non-comoving multiple 
systems 

HD 3003 — /3 3 Tucanse (HD 3003) was previously known as 
a multiple star (CCDM J00327-6302AB). Were this system co- 
moving, we should be able to retrieve both components, sepa- 
rated by 674 in 1925 (as tabulated in the CCDM). Only one star 
is visible in the 2005 NaCo image, hence this system is only a 
visual assocation. 



HD 50445 — The star HD 50445 is not catalogued as being part 
of a multiple apparent or physical system. However, it lies at a 
relatively low Galactic declination (b = -15.6°) and at least 5 
point sources have been detected in 2005-November in the 27"- 
wide field of NaCo S27 camera, using the coronagraphic mask. 
All 5 objects were observed again in 2009-April; a software is- 
sue during the observations unfortunately prevented recording 
a direct image of the primary star at this epoch, hence the as- 
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trometric precision is lower for the 2009 observations (typically 
1 pixel, i.e., 30 mas). We are nevertheless confident that all can- 
didate companions are background contaminants. Note that CC2 
can be seen in the infrared Digitalized Sky Survey 2 image taken 
on 1985-December atp ~ 14" and 6 ~ 155° (the apparent mo- 
tion of the background contaminants on the sky between 1985 
and 2005 is /x a ~ +078 and \i s ~ +172). 



(a) Confirmed companions 




HD 1 77756— A Aquilae (HD 177756) is one of the B stars in- 
cluded in the Hubrig et al. (2001) survey. These authors did not 
report any companions in the K band with a magnitude differ- 
ence up to AK ~ 9. The observed CC1 has a magnitude differ- 
ence of AK - 4. It should have been seen in the 1999 ADONIS 
image, as this background contaminant was closer to the star at 
this epoch, unless the field of view of 24" x 24" reported by 
Hubrig et al. is not squarred (depending on the dithering pattern 
used by these authors). 

HD 224392 — Members of the Tucana association, including 
the Vega-like star r\ Tucanse (HD 224392; see Mannings & 
Barlow 1998 for the infrared excess), were surveyed in the in- 
frared by Neuhiiuser et al. (2003) with the SHARP -I infrared 
imager on the New Technology Telescope (NTT) at La Silla. 
These authors report detection limits in the K band of AK - 8.3 
at > 100 au from the star. However, they could not have de- 
tected CC1 for it was out of SHARP-I 13" x 13" field-of-view. 
Zuckerman, Song & Webb (2001) note that as part of the IRAS 
Faint Source Catalogue (Vizier online data catalogue II/156A; 
Moshir 1989), r\ Tuc might also be a multiple star, yet in this 
case we do not resolve it. Besides, being a HARPS constant 
(Lagrange et al. 2009a) makes it unlikely to be a spectroscopic 
binary. We do detect two companion candidates far from the star, 
which are not comoving with it. 



5. Discussion 

5.1. Multiplicity of early-type stars in the sample 

We found that 6 stars out of 38 in the sample have confirmed 
low-mass companions, yielding an observed multiplicity frac- 
tion of 16% for our sample. This number can be seen as a lower 
limit value for the sample (i) because we have imaged compan- 
ion candidates around 7 other stars in the sample, without having 
the possibility to determine whether they are true companions 
or background objects, (ii) because we could have missed close 
binaries below our detection limits, and/or (iii) because some 
wide-separation binary components could have been too close 
to the line of sight to be resolved at the time of observations. It is 
also certainly a lower-limit value regarding the global multiplic- 
ity fraction of early-type stars, since our survey is biased against 
spectral binaries or visual binaries with companions at < 5", as 
described in Sect. 2.1. 

In the case all the undefined candidates are bound to their 
primary stars, the multiplicity fraction of sample stars could be 
as high as 34%. More specifically, the multiplicity fraction for 
F stars in the sample is 5/16 ~ 31%; it is 1/22 w 5% for 
A+B stars. If we assume that unconfirmed companions are all 
physically bound companions, then these fractions would rise 
to 50% and 23% for F and A+B stars, respectively. The multi- 
plicity fraction could be even higher than these figures due to 
narrow-separation binaries and orbital projection effects. 

Considering confirmed companions only, it would seem sur- 
prising that the multiplicity fraction of F stars is 6 times larger 



(b) Confirmed + unconfirmed companions 
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Fig. 12. Number of detected companions as a function of the projected 
separation to the primary, per bin of 100 au. Companions around F stars 
(red stripes) and A+B stars (green stripes) are differenciated. The black 
dotted line represents the total number of companions, (a) Confirmed 
companions only, (b) Confirmed and unconfirmed companions. 



than that of A stars, whereas it is usually expected that the mul- 
tiplicity fraction increase with the stellar mass. Actually, this is a 
bias that can be attributed to the fact that we are probing different 
physical separations for different stellar types. In order to have 
about the same number of F and A stars included in this volume- 
limited survey, we have to go farther to include as many A stars 
as F stars: the median distance of sample F stars is 27.0 pc while 
it is 54.9 pc for sample A and B stars (see Fig. 3). Thus, we are 
probing smaller physical separations for F stars than for A and 
B stars, as shown in Fig. 7b. Five out of 6 confirmed compan- 
ions detected next to F stars are located between 10 and 40 au, 
at projected separations not or hardly reached for the A stars in 
the sample. On the contrary, a confirmed companion located at 
> 100 au was found for each type of stars. 

Figure 12 shows the number of confirmed companions 
(Fig. 12a) or confirmed and unconfirmed companions (Fig. 12b) 
per bin of 100 au from their primary stars. The bias is clearly 
apparent in Fig. 12a if we assume that our detection efficiency 
at < 100 AU is lower for A stars than for F stars. If we further 
assume that all unconfirmed companions are physically bound 
to their primary stars, then Fig. 12b shows that we have roughly 
the same number of objects between 100 and 500 au, while com- 
panion candidates at projected distances greater than 500 au are 
only found around A+B stars. It is likely that at least some un- 
confirmed objects detected around A stars are real companions, 
and contribute with the observational bias to fill the gap between 
the multiplicity fractions of A and F stars in the sample. 

5.2. Impact of imaged companions on the velocimetric 
survey 

The objects from the imaging survey presented in this work have 
been - or are being - monitored spectroscopically, seeking for 
radial velocity variations that would reveal the presence of close- 
in companions (see Lagrange et al. 2009). Two questions arise 
when companions are imaged around stars included in the ve- 
locimetric survey. 

First, is the spectroscopic signal contaminated by the pres- 
ence of a companion that would appear blended with the primary 
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star in the spectrograph fiber? All background objects detected 
in the imaging survey are far enough from their stars and do not 
enter the fiber (© ~ 2") used in the velocimetric survey. Thus, 
they do not pollute the radial velocity signal. On the contrary, the 
companion candidate to HD 43940 (ambiguous case) probably 
contribute to the velocimetric signal. Six of the seven comov- 
ing companion detected, HD 14943B (p = 273), HD 49095BC 
(p = 173), l Crt B (p = 079), 26 Oph B (p = 073), and o Gru B 
(p = 074), are close enough to their stars to directly contribute to 
the spectroscopic signal. Note that o Gru B probably marginally 
pollutes the velocimetric signal depending on seeing conditions. 
However, the contribution in terms of flux is generally limited 
given the flux ratios between the stars and the companions. 

On the other hand, are the radial velocity variations a result 
from the gravitational perturbations of the imaged companions? 
Indeed, all the 6 stars for which we imaged and confirmed a com- 
panion are reported bearing velocimetric variations in Lagrange 
et al. (2009). Radial velocity curves (Lagrange et al. 2009) in- 
dicate that two of them (HD 41742 and o Gru) belong to mul- 
tiple systems; this includes a spectroscopic binary surrounded 
by a wider-separation companion, detected only with adaptive 
optics. 6 The imaged companions are not responsible in these 
cases for the radial velocity variations. Meanwhile, the con- 
firmed companions around the other 4 stars partly or totally con- 
tribute to the radial velocity variations observed by Lagrange et 
al. (2009) as well as in new unpublished data (Lagrange et al. 
in preparation). In the ambiguous cases of HD 43940 and rj Hor, 
the companion candidates could contribute to the observed radial 
velocity variations, if they are bound objects. 

In any case, this shows that the radial velocity curves of at 
least 6 out of 35 stars followed spectroscopically (17%) may be 
affected by stellar multiplicity. 

6. Conclusions 

We have presented a deep imaging near-infrared survey of 38 
austral early-type main sequence stars, searching for low-mass 
companions at intermediate distances. The overall detection per- 
formances obtained for this survey, mainly performed in the K 
band at 2.2 //m, allowed us to probe the stellar environments 
at separations > 0/2 with contrasts < 0.05. Farther from the 
stars (> 5"), we achieved limiting contrasts between 10~ 5 and 
10~ 6 using coronagraphy. Although known binary stars were ex- 
cluded from the survey, we detected 41 companion candidates 
around 23 stars. Using multi-epoch observations, we were able 
to determine that 8 candidates are actually comoving compan- 
ions around 7 stars of the sample, while the other candidates are 
mainly background contaminants. The comoving companions 
are low-mass stars with masses ranging from 0.13 to 0.8 M Q . 
A comoving binary system consisting of 2 M/L dwarfs imaged 
around HD 49095, with an almost fully resolved orbit, will be 
of particular interest for the calibration of low-mass star evo- 
lution models. The comoving companions are detected at pro- 
jected separations ranging from 1 1.3 to 159 au. 

At least 16% of the sample stars, with spectral types rang- 
ing from F7 to B9, therefore host a low-mass star companion. 
This figure is likely a lower limit to the real multiplicity frac- 
tion of early-type stars within this separation range. In addition, 
all the detected and confirmed companions have projected sepa- 
rations below 200 au: this range matches the physical projected 
separations of giant planets recently imaged in discs surrounding 



6 In the case of HD 41742, there is a 4th companion lying out of the 
field of view of this survey. 



A-type stars Fomalhaut and HR 8799 (Kalas et al. 2008; Marois 
et al. 2008). It is certainly tantalizing to interpret the emerging 
population of such massive planets as the low-mass tail of the bi- 
nary distribution of early-type stars. This would imply that such 
planets are formed by gravitational instability rather than core 
accretion like closer-in planets. However, further evidence from 
larger, less-biased surveys, are still required to assess this ques- 
tion in a statistical and quantitative fashion. 

Most stars of our the sample are also included in spectro- 
scopic surveys (Lagrange et al. 2009) seeking for radial velocity 
variations. Crossing the results of both surveys show that at least 
17% of the radial velocity curves may be affected by companions 
detected with adaptive optics. 

We finally emphasize that tight stellar systems with a bright 
early-type primary component, such as the ones presented in 
this study, will be excellent calibrators for next-generation 
planet imagers such as VLT/SPHERE or interferometers like 
VLT/Gravity, designed for high-precision astrometry within 
small field-of- views. 
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Fig. 4. Examples of the status determination of companion candidates around HD 14943 (left) and rj Tuc (HD 224392; right), from astrometric 
measurements at two different epochs. The positions on the sky are reported in the upper panel: the candidates were first imaged on 2005-11- 
07 (filled black circles), then on 2008-08-20 (filled red and blue circles; there are two measurements in the case of HD 14943). The positions 
predicted for the candidates on 2008-08-20 if they were background objects from the field are indicated by the empty red and blue circles, and 
their trajectories between the two epochs figured by the black curves, which take into account the stellar proper motions and the parallactic 
motions. The middle and lower panels show the evolutions of the projected separations and the position angles, respectively, with time. One-sigma 
error bars are represented. A zero or near- zero evolution of the separation and position angle with time is indicative of a physical association. The 
candidate at the left is bound to its star with indicative probabilities P^kg = 0% and P cmv = 96.7% (following the calculations of Sect. 4.2), while 



the candidate at the right is a background contaminant, with P bkg = 12.0% and P a 



0%. 
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Fig. 5. Example of 6-<r detection limits reached for a K = 4.7 mag A0V star at 52.7 pc (HD 2834) using direct imaging (left panel) and coronagra- 
phy (middle panel). Detection limits for the direct and coronagraphic image are given in 2 x 2 and 6x6 arcsec 2 fields around the star, respectively. 
The hatched area on the coronagraphic image represents the occulting mask with a diameter of 077. Classic 1-dimensional 6-cr detection limits 
(right panel) are extracted from the direct image detection limit map (plain line) and the coronagraphy detection limit map (dotted line) as ex- 
plained in Sect. 4.3. Mass limits (red lines) are obtained by interpolating Baraffe et al.'s (1998) evolutionary model for a given stellar age (here, 
700 Myr). 
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Fig. 6. Same as Fig. 5 f or a K = 4.3 mag F7V star at 24.6 pc, and an age estimated to - 4. 1 Gyr (HD 9 1 
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(a) HD 14943 (b) HD 41742 
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Fig. 9. Stars with companion candidates confirmed through multi-epoch observations. North is up and east is to the left. Counts on the detector are 
displayed with a logarithmic scale. 
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Star name Proper motions Spectral Distance Apparent magnitudes Age Multiplicity RV oo 

a S type V J H K ref. code ref. code ref. 

(mas yr~ ' ) (pc) (mag) (mag) (mag) (mag) (Myr) 

Table 1. Star sample for the southern survey. Infrared J, H, and K magnitudes are extracted from 2MASS (Skrutskie et al. 2006). The 'RV column 
indicates the radial velocimetry status (C: constant, V: variable, with possible source of variations being SB: spectral binary, A: magnetic activity, 
Pu: pulsations, PI: planets, D: drift) as determined by Lagrange et al. 2009 with HARPS or as relying on SOPHIE data (Desort et al., personal 
communication). The column 'Multiplicity' indicates apparent (or comoving) binary (B), ternary (T), or more-component (M) systems, whether 
it is new (N) or known before this survey (K). References are given below the table. 

References: (1) Su et al. 2006; (2) HARPS survey, Lagrange et al. 2009; (3) HARPS survey, unpublished data; (4) CCDM II, Vizier Online Data 
Catalogue 1/274, Dommanget & Nys 2002; (5) Fernandez, Figueras & Torra 2008; (6) Geneva-Copenhagen Survey of Solar neighbourhood III, 
Vizier Online Data Catalogue V/130, Holmberg, Nordstrom & Andersen 2009; (7) Multiple Star Catalog, Tokovinin 1997; (8) Proper motion 
derivatives of binaries, Makarov & Kaplan 2005; (9) SOPHIE survey, unpublished data; (10) S B o, Pourbaix et al. 2004; (11) Desort et al. 2008; 
(12) Goldin & Makarov 2007; (13) Eggen 1998; (14) proper motion binary, Frankowski et al. 2007; (15) Rieke et al. 2005; (16) The Washington 
Visual Double Star Catalogue, Mason et al. 2001; (17) King et al. 2003; (18) member of the Tucana association, Zuckermann, Song & Webb 2001; 

(19) HD 4293 is classified as a F0V but has colours corresponding to a F0III giant (Hauck 1986) with typical age between 10 and 100 Myr (Eggen td 

1991a); (20) Morales et al. 2009; (21) Plavchan et al. 2009; (22) Gray et al. 2006; (23) Makarov 2007 suggested that HD 19545 and BO Microscopii gf 

('Speedy Mic') were formerly part of the same system within the Tucana-Horologium association; (24) Grosb0l 1978; (25) member of the IC 2391 H 

supercluster (Eggen 1991b); (26) a field B star, Huang & Gies 2008 o' 
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Telescope/Instrument 
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VLT/NaCo 
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VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 
CFHT/PUEO 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 
CFHT/PUEO 
CFHT/PUEO 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 
CFHT/PUEO 

VLT/NaCo 
CFHT/PUEO 

VLT/NaCo 
CFHT/PUEO 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 
CFHT/PUEO 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo 

VLT/NaCo __^__ 

Table 3. Observing setups, dates, and number of companion candidates observed ('M' stands lor many). For each observing date, we precise the 
observing mode (D: direct imaging, C: coronagraphic imaging, S: saturated direct imaging), the filter(s) (NaCo filters: J, H, K s , NB2.12; PUEO 
filters: K', Fen, Bry), and camera(s) used (NaCo cameras: S27, S13; PUEO camera: KIR). The last column 'Companionship' indicates the 
presence of companion candidate(s) (CC) or bound companion(s) (BC) when this has been established at several epochs. The epochs II, V, VI, 
and VII correspond to ESO observing programmes 076.C-0270(A), 081.C-0653(A), 083.C-0151(B), and 083.C-0151(A), respectively. Epochs I, 
III, and IV correspond to CFHT programmes . . . 
*No second-epoch observation was recorded. 
**No coronagraphic image was recorded. 
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Star 



CC# 



Date 



Filter Camera Mode 



DIT X NDIT X NEXP 
(s) 



Projected separation 
P(") 



Position angle 

en 



Am ± 0.4 
(mag) 



Status 



10. x 3 x 20 
10. x 3 x 20 

40. x 1 x 15 
2.0x3x 10 

0.5 x 1 x 25 
0.5x20x 12 
0.3x20x 12 

2.0x5x 10 
0.8 x 5 x 60 
0.8 x 5 x 20 
2.0x5x 10 
0.8 x 5 x 60 

0.5x60x 10 
0.5 x 60 x 10 
0.5x60x 10 

1.0 x 10 x 10 
1.5x20x8 

0.8 x 30 x 20 
0.5x80x 10 

3.0x9x20 
5.0x6x20 
2.0 x 15 x 10 

0.0 x 1 x 30 
1.5x20x 10 

4.0x8x 10 
2.5 x6x 10 

0.8 x 5 x 10 
1.0 x 10 x 10 
4.0 x 3 x 10 

0.0 x 3 x 20 
0.0 X 3 X 20 
4.0x8x 10 
4.0x8x 10 

Table 5. Multi-epoch measurement results for 41 companion candidates detected around 23 stars of the sample. The imaging mode used for each 

detection is indicated: (D) direct imaging, (S) saturated imaging (only on PUEO), (C) coronagraphic imaging (only on NaCO). The companion 

candidate status are defined in Sect. 4.2. 

Notes for some targets: 

HD 32743 — p, 6, and Am of CC#2, 3, & 4 have been determined with respect to the properties of CC#1, which is seen both in direct and 

coronagraphic images. 

HD 153363 — Only the two brightest companion candidates are reported. Many other point sources (~ 56) are seen using the coronagraph; 

however, they are all likely background contaminants as the star is close to the Galactic plane (b = +10.6°). Besides, we miss a second epoch 

observation with the coronagraphic mask. 

HD 216385 — CC#1 is not seen in the second epoch observation, while CC#2 is too weak to have been detected during the first PUEO epoch. It 

was, however, out of the field of a third NaCo epoch not reported in this Table, and is thus a likely background contaminant. 
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